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, A generalized scattering amplitude where momenta of incoming-particles and outgoing- 

particles as well as positions of incoming-particles and outgoing-particles are specified is 
formulated. Idealistic beams and idealistic measuring instruments where momenta and posi- 
| tions satisfy minimum uncertainty are studied with a use of minimum wave packets, coherent 

states. In the present work, we show general features of the generalized scattering ampli- 
tudes based on cjf> 4 theory. We give a proof of completeness of many body states, asymptotic 
behaviors in the large distance region, and factorization of the amplitudes. Despite of the 
non- orthogonal properties of wave packets, we found that the probability interpretation is 
verified. A differential probability depends upon the wave packet size but a total probability 
that is integrated in the final states is independent from the size of final state wave packet 
and becomes universal. Few body amplitudes are studied as examples. 
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§1. Introduction 
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We study scattering processes in which both of momenta and positions are si- 
multaneously measured. Observations of both variables are made in recent neutrino 
experiments where a distance between a source region and a detector is fixed in 
a macroscopic length and the momenta of the particles are measured. Scattering 
amplitudes of both variables are formulated in the present paper. 
Ph, In quantum mechanics, precise values of momentum and coordinate of a particle 

Q_ij can not be measured simultaneously from Heisenberg uncertainty relation. So in 

scattering processes, one of the variables is selected and a dependence of the scatter- 
ing amplitude on this variable is studied. Normally momentum variable is selected 
and the momentum dependence of the transition probability is studied and is com- 
^ ■ pared with a theoretical calculation. Since momenta are commuting variables, it is 

possible to determine them precisely. In real experiments, however, the exact value 
of momentum is difficult to measure and the momentum is measured with finite un- 
certainty due to a finite resolution of a detector. When the momentum is measured 
with uncertainties, simultaneous measurements of positions and momenta become 
possible in experiments. The positions are measured also within finite uncertainties 
due to spatial resolution of detectors. When momenta and positions are measured 
with finite uncertainties, total number of information could be equivalent with the 
standard case where the exact values of momenta are measured. 

Since a wave function with a definite momentum is a plane wave and is invariant 
under any translations, any particular position is undefined. But a wave function 
with a finite uncertainty of the momentum can have a finite spatial extension and 
is localized around certain position, so in this case the spatial position is defined. 
So it is necessary to introduce finite uncertainty of momenta in order to define the 
positions. As a price of a finite uncertainty on the momentum, it became possible 
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to introduce a position where the particle is measured or is produced. The position 
has also finite uncertainty. In this work, we introduce such amplitude that is defined 
with finite uncertainties of momenta and positions. 

The wave functions with finite spatial extensions are necessary actually for 
asymptotic conditions of scattering processes to be satisfied 1 ) . 2 ) However in many 
situations of high energy physics, effects of finite wave functions have been ignorable 
and it has been sufficient to use propagators of momentum variables with ie pre- 
scription. We clarify these points and we study situations where the effects of finite 
wave functions are important. To investigate these problems we define the gener- 
alized S-matrix of both variables using mixed representation and we show several 
features of the scattering amplitudes in which the momentum and the positions are 
measured. Some implications and applications are also studied. 

We formulate a position-dependent and momentum-dependent S-matrix in an 
idealistic limit allowed from Heisenberg uncertainty relation. In our formalism, each 
momentum and position of the initial particle and final particle satisfy minimum 
uncertainty relation. This corresponds to a scattering process where coordinates as 
well as momenta of the beam satisfy the minimum uncertainty relation and the mea- 
surement are made with minimum uncertainty allowed from Heisenberg uncertainty 
relation. 

In real experiments the uncertainty of the momentum and coordinate may be 
larger than those of the present work. But our idealistic scattering matrix should be 
realized in a suitable method and should give new insights in quantum mechanical 
scattering processes. Also extensions to non-minimum wave packets is straightfor- 
ward. 

Uncertainties of momenta and positions are actually determined in experiments 
by resolutions of beam sources and of detectors. The resolutions depend on each ac- 
celerator and detector. So we leave the uncertainties unfixed and study the scattering 
matrix of this situation. 

Because the momentum is conjugate to the coordinate, they satisfy the commu- 
tation relation, 

[x,p]=i, (1-1) 

in a unit h /2ir = 1 and a momentum resolution and a spatial resolution of a detector 
satisfy Heisenberg's uncertainty relation, 

S p x Sx > (1-2) 

A product of uncertainties becomes minimum in coherent states. The coherent 
state of a variable, x, 

(x\P ,X ) = jVV^-Xo)-^-^) 2 ^ 
Nf = (trt)-5 (14) 

has expectation values 



<M> = *o 



(1-5) 
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(|x 2 |) = X 2 + ia 
(\p\)=Po 
(\P 2 \) = Po 2 + 



(1-6) 
(1-7) 
(1-8) 



The product between the variances of the momentum and coordinate, 



(fe) 2 x (<5p) 2 = \ 

(5x) 2 = (\x 2 \) - ((\x\)) 2 
(5p) 2 = (\P 2 \) ~ ((\P\)) 2 



(1-9) 

(MO) 
(Ml) 



is independent from a and satisfies the minimum uncertainty condition. The uncer- 
tainty becomes minimum with the coherent state. 

The coherent state satisfies also completeness condition, 



We use these wave functions for expressing in-states and out-states. In the 
former, the state corresponds to the idealistic beam and in the latter the state 
corresponds to the idealistic measuring apparatus. A transition element between the 
states of idealistic beams and the states of the idealistic measurement is computed 
from these matrix elements. Using them, we define the idealistic scattering matrix 
in mixed representation and study its properties and applications. 

Our formalism will be applied to long line experiments such as, solar neutrino 
experiments, atmospheric neutrino experiments, long base line neutrino experiments, 
reactor neutrino experiments, and others where systems have huge scales and posi- 
tions of detectors play important roles, SM)^ 8 ) 9 ) 10 ). 13 ) i n addition to neutrino, scat- 
tering of other weakly interacting particles where position dependence in addition to 
momentum dependence are measured and give important physical informations will 
be studied. 

The paper is organized in the following manner. In Section 2, mathematics 
on wave packets are given. Wave packet wave functions are explicitly given and 
the overlap integrals and time dependent behaviors are analyzed. New uncertainty 
relations between the velocity of expansions in asymptotic region and the initial 
sizes are obtained. Although some of the materials in this section may be known to 
experts, they are necessary and useful for later arguments. In Section 3, generalized 
scattering amplitude is defined and general properties are studied. We find suitable 
integration measures for both variables in which the probability interpretation is 
verified. A differential probability depends upon the wave packet size but a total 
probability that is integrated in the final states is independent from the size of final 
state wave packet and becomes universal. Explicit examples are given in Section 4 
and summary is given in Section 5. 




J 27T 

= S(x - y). 



(M2) 
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§2. Wave packets 

2.1. Mathematical preliminaries 

2.1.1. Complete sets of minimum wave packets :time independent wave packets 

Let momentum eigenstates and position eigenstates be \p) and \x). Transforma- 
tion from a momentum eigenstate to a position eigenstate is made from 

(x\p) = (27r)-§e** (2-1) 
<£|p>* = (p\x) = (x\-p). (2-2) 

Both sets of functions are complete sets and satisfy 

dx(pi\x){x\p 2 ) = 6(pi - p 2 ), (2-3) 

dp{xi\$(p\x 2 ) = 6(xi - x 2 ), (2-4) 

and 

' dx\x)(x\ = l, (2-5) 

dp\p)(p\ = 1. (2-6) 

The normalized coherent states in three spatial dimensions 1 ^ are defined as 

(x\P , X ) = N 3 e iPo ^-^-^ s -^ 2 (2-7) 

(x\P ,X )* = (P ,X \x) = (x\-P ,X ). (2-8) 

Due to the Gaussian factor the wave function is localized around the center position 
xq with a width \f 7 lo and is an approximate eigenfunction of the momentum operator. 
A constant N% is determined from the normalization condition 



/' 



/ 



d 3 x\(x\P ,X }\' 



No 



J d 3 xe~^-^ 2 (2-9) 



= iVf (7TCj)5 
= 1. 



The same states are expressed in the momentum representation as 
(p\P ,X ) = J dx(p\x)(x\P ,X ) 

= N 3 a 3 / 2 e- l P-^-^ p °-?> 2 , (2-10) 
(p\P ,X )* = (P ,X \P) = (p\-P ,X ) (2-11) 
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and satisfy the normalization condition in the momentum representation. 

The set of functions for an arbitrary value of a satisfy the completeness condition 
in the coordinates representation, 



/ 



^^(x 1 \P ,X )(P ,X \x 2 ) = 5(x 1 -x 2 ), (2-12) 
and the completeness condition in the momentum representation, 

/ ^^(Pi\PoJo)(Po,Xo\P2) = m-P2)- (2-13) 
The completeness condition is simply expressed as 



/ 



^^\P ,X )(P ,X \ = 1. (2-14) 

Overlap between two coherent states are computed from the above definition 
(P 1 ,X 1 \P 2 ,X 2 ) 
= J d 3 x(P 1 ,X 1 \x)(x\P 2 ,X 2 ) 

= N$ y*d 3 xe ii£ ^ l3:+i( ^ Xl - P2X2) -^( ( ^ l)2+( ^ 2)2 ) (2-15) 

= e -^(^i-^2) 2 -f(Pi-P 2 ) 2 + f(Pi+P 2 )(X 1 -X 2 )_ 

Obviously matrix elements do not vanish and states are not orthogonal for dif- 
ferent values of the momenta and coordinates. Despite of the nonorthogonality they 
satisfy 

dP dX . 



I 



(27 _ ):{ \P ,X )(P ,X \P,X) = \P,X). (2-16) 

Hence (P\, Xi\P 2 , X 2 ) plays a role of the Dirac delta function. 
2.1.2. Wave packets defined at arbitrary time 

Wave packets defined at a certain time To is constructed when one particle energy 
is known. Let E(p) stands one particle energy of the momentum p, 

E{p) = (f + m 2 ) 1 / 2 (2-17) 

in the unit with c = 1, then the wave packet defined at To is, 

(p\P ,X ,T ) = (p\P ,X )e^ T °. (2-18) 
This set of functions for a given time To satisfy the completeness condition, 
dP dX 



I 



(2tt)S 



■(pi|To,^o,To)(To,Xo,To|p 2 ) = (n\p 2 ). (2-19) 
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The representation of these wave packets in the coordinates space is obtained by the 
Fourier transformation, 

(x\P ,X ,T ) = J d 3 p(x\p)(p\Po,X )e^ T °. (2-20) 

It is easy to verify the completeness condition by combining Ea. (|2-14[) and (|2-18|) 
with Ea. (t2~20|) . 

/ ^^ (^\Po,X ,To)(P ,Xo,To\x 2 ) = (xi|z 2 )- (2-21) 
The completeness condition is simply written as 

/ ^{itf ° \Po,Xo,T ){P ,X ,To\ = 1. (2-22) 
The matrix elements of the wave packets at different times are computed as, 

(P 1 ,X 1 ,T 1 \P 2 ,X 2 ,T 2 ) = j d 3 p{P 1 ,X 1 ,T 1 \p)(p\P2,X2,T 2 ) (2-23) 

and satisfy 

dP dX ,~ 



1 



^_ y> f P ,X ,T )(P ,Xo,T \P 1 ,X 1 ,T 1 ) = \Pi,X 1 ,T 1 ), (2-24) 

{PuXx^P^X^Tx) = (P^jt^Xt). (2-25) 

Explicit forms of these matrix elements are given later. We use these wave functions 
for expanding the field operator. 

2.2. Matrix elements 

2.2.1. Time dependent transformation function 

We calculate the matrix elements of the mixed states appeared in the previous 
section. It is convenient to define time dependent transformation functions, 

(t,p\P ,X ,T ) 

= e E ¥ l ( t - T °Hp\P ,X ) (2-26) 



and 



(t,x\P ,X ,T ) 

= J dp(x\p)(t,p\P ,X ,T ) (2-27) 

2ir J 
2vr J 
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The absolute value of the integrand becomes maximum at Pq but the phase becomes 
large in large t — Tq region or large x — Xq region. So we integrate on p in two 
regions separately. In the small t — Tq region we use the approximation of the 
integrand around the Po an d in the large t — Tq region we use the approximation of 
the integrand around the stationary momentum. 

(A) Small T — Tq case:translational motion. 

In the small t — Tq region, the integral is written and computed around Po as 
(t,x\P ,X ,T ) 

= jV 3 (-^) 3 / 2 f ^ e -»(^(Po)+(p-po)-w)(t-r )+j(po+(p-po))<?-Xo)-|(p-Po) 2 
2ir J 

= Ne**, (2-28) 
N = N3e -^(S-x -vo(t-T )) 2 ^ ( 2 -29) 
e i<t> = e -iE(Po)(t-T )+iPoiS-Xo) ^ (2-30) 

^ = Wi Emp=P ~°' (2 ' 31) 

The wave packet keeps its shape and moves with a constant velocity vq. The center 
of wave packet is Xq at t = Tq and is Xq + vo(t — Tq) at a time t. 

(B) Large t — Tq: expanding wave packet. 

In the large t — Tq region the momentum integration applied in the previous 
method is not a good approximation any more because the phase oscillates rapidly 
in this region. The phase of the integrand becomes stationary at Px which satisfies 

^- [-iE(p)(t - T ) + if. (x-Xq) - G -{p- P ) 2 ) = 0. (2-32) 
The solution is obtained by expanding the momentum in -j-jjr and is given by, 

Px = 4 0) +4 1) + 4 2) , (2-33) 

P| 0) = m 1 (x-Xq), (2-34) 

^(t-T )2-(x-X ) 2 

Px = i-^frvE{P x ){P x - P ), (2-35) 

t — lQ 

4 2) = O ((t-T )- 2 ). (2-36) 
So the exponent of the integrand is expanded around Px and we have 

(t,x\P ,X ,T ) (2-37) 

= N e - iE tP){t-To)+iV-{x-Xo) 

N = N 3 (-^— )V^-±— )e -\^ ) -^\ (2-38) 



K(Px-Pq) 2 {2i-^ (2-39) 
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( = P%1, (2-40) 

(t-r )' 

in a large t — Tq where the wave packet parameters are given by 

lm 2 \t-T \ 

1L = — — -, 2-41 

2 E {p x y 1 ] 

lit — XI 

7T = - J — ^- L - (2-42) 

The longitudinal component and transverse component of a momentum q are defined 

as 

The phase factor in the Eq. ()2-37j) is written in leading order of \t — Tq\ as 

e -i{t-T )E(P x )+iP x -(x-Xo)) 

= e -i^(t-T )^-{s-x Y (2-45) 

(t-To) 



= e E< - p xi 

This phase factor becomes very small if the mass is very small and vanishes in the 
massless case. The fact that the phase becomes small in the high energy region or in 
the small mass region is a characteristic property of relativistic invariant theory. 12 ) 
The absolute magnitude of N becomes maximum when the momentum Px 
agrees to Pq. This is realized at a particular position of f , xq, which satisfies 

P = m — 1 (xp - X ). (2-46) 

y (* — ^o) 2 — (%o — Xq) 2 

The solution is 

x = X + (t-T )-^- (2-47) 
ti{Po) 

where xq is the center of wave packet at a time t. 

We write the difference of the momenta, Px — Po, m the Gaussian exponent in 
Ea. (|2-37|) using x and xq and a unit vector n\ in the pi direction as, 

Px ~ Po = ^-{(x - x )t + - fo)L|ni}. (2-48) 

We substitute this expression into Eq. (j2-37|) and we have 

exp(~a(P x - Pof) (2-49) 

1 (E(P ))\ ^ 2 1 (E(P )f _ 2 

= eM -2 a JT-T T {x - Xo)T -2 a m*(t-T r {x - Xo)L } - 
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The normalization factor in the small \t — Tq\ region has also 

exp(-^(x-x ) 2 ). (2-50) 
Hence the size of the wave packet in the longitudinal direction is given as 

6x L = JIH^M + ^ ( 2 .5l) 
Vex E(P f 

and in the transverse direction is given as 

2\t-T \ 



Sx T = J-i + V2a. (2-52) 

V o E(P ) 

A wave packet expansion is characterized by its velocity. The velocity of expansion 
in the transverse direction, vt, is determined by the momentum variance as 

v T = \-—^. (2-53) 
V^£(Po) 

and that in the longitudinal direction,^, is determined as 

v L = (2-54) 
V a {E {P )f 

The velocity of expansion satisfy uncertainty relations 

v T 5x(t = 0)E(P ) = 1 (2-55) 
v L 5x(t = 0)E(P ) = (-^-) 2 , (2-56) 

where 5x(t = 0) is the spatial extension of wave packet at t = 0. The vl is given 

2 

by multiplying ^ to the vt- This ratio becomes one in the non-relativistic 

energy region where E(po) is nearly equal to m. Consequently wave packets expand 
symmetrically in the non-relativistic region. In the relativistic region where E(po) 
is much larger than m, the ratio between both values becomes very small and wave 
packets expand an-symmetrically. The shape becomes a circular thin disk after 
certain time. The wave packet size is given as a function of the propagation time, t, 
for various values of the initial wave packet size in Fig(l). The 5xt becomes huge 
size in 500 sec, the period between the sun and the earth. 

Due to the expansion of wave packet, the normalization of the wave function is 
inversely proportional to 7t\/7l in the asymptotic region. 

2.2.2. Two body matrix elements 

Various matrix elements of wave packets are studied in this section. 
1. Matrix elements of wave packets defined at equal time 
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10" 



10° 



10 u 



10" 



10"' 



1 1 1 1 1 1 1 1 , I J. J. J- J- 

+++++++++ +++ + ++ + ++++++++++++++++++++- u+++ 

+ + + + + + vv vxxxxXXXXXXXXXXXXXXXXXXXXXXXXX 



X xxxxxxxxxxxxxxxxxxxxxx 
, ^ X * ^ ^ 



x T , 0=1.0*10" 
x T , o=1.0*10" 1 °m 

x T , 0=1.0*10", m 
x Ll o=1.0*10 " m 
x L , a=1.0*10" gm 

x L , 0=1.0*10 m 



□□□□□□□□□□□□□□□□°° DnDnnnnDnDDDDDDnDDDDDDI1 
□ DDDDD 

giB8oS5oOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 
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Fig. 1. Time dependence of wave packet size in the physical unit is given. St is the size in the 
transverse direction given in Eq. l2-52H and Sl is the size in the longitudinal direction given in 
Eq.lJjSnj • mc 2 = 10" 2 eV is assumed. 



The overlap between two wave packets is given as, 

(P 1 ,X 1 \P 2 ,X 2 ) (2-57) 
= e -2^i-^) 2 -f(Pi-P 2 ) 2 e i(Pi+P 2 )(Xi-X 2 ) _ 

Thus the overlap decreases fast as the distance between two positions in the real space 
and in the momentum become large. Matrix elements for the same coordinates is 

(Px,X\P2,X) = e--^- p rf (2-58) 

and for the same momenta is 

(PJilPjz) = e -^(^-^) 2 +iPi-(x 1 -x 2 )_ (2 . 59) 

From these matrix elements if a particle is prepared at (P',X') probabilities of 
finding a particle in a region of the momentum and coordinate between (P, X) and 
(P + dP,X + dX) is given by 

P = ^L e -^(x-xr-^P-P^ (2 . 60) 

(27r) ,i 

Heisenberg uncertainty relation between the variance in the coordinates and the 
variance in the momenta is satisfied. It is important to notice that the momentum 
and the coordinate is measured simultaneously. 
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2. Time dependent matrix elements 
The time dependent matrix element is computed by inserting the complete set of 
momentum eigenstates, 

(Pl,Xi,Ti\P 2 ,X 2 ,T 2 ) 



I 



d 3 p(P 1 ,X 1 \p^e ^{p\P 2 ,X 2 ) (2-61) 

2(J 3 e -f (A-P2) 2 y d 3 pe EmT ? " T2) +»p-(X 1 -X 2 )-a(p-P ) 2 ; 



P = 1(^+^2). 

In the small \T\ — T 2 \ region this matrix element is given as, 

(Pi,^i,Ti|A,X 2 ,T 2 ) (2-62) 
= iY e - iS ( /3 o)(Ti-T2)+iPo-(X 1 -X 2 ) ) 

where the velocity is given by 

«_«M|,_Ju (2 . 64 ) 

This matrix element shows that the shape of the wave packet is preserved but the 
center position is moving with the velocity v in the small \T\ — T 2 \ region. 

The matrix element in a large \T\ — T 2 \ region is obtained by the stationary 
phase approximation in p integration. The stationary point Px is obtained from the 
stationarity condition, 

— (if- (Xi - X 2 ) - iE($){Ti ~ T 2 ) - P ) 2 ) = 0, (2-65) 

OPi 

and is given as, 

P x = <X 1 - X 2 ) ^ ^ + i 5P X , (2-66) 

5P X = 2aE(P x ){P x - Po). (2-67) 
Using this momentum, we have 

{Pi, X\,Ti\P 2 , X 2 ,T 2 ) 

= Nexp(i(j)), (2-68) 

N = (1 + ^-JV(1 + ^-T 2 W /2 
2oE{P x )' 2aE{P x y 

xe -f{(A-P 2 ) 2 +4(P x -P 0) 2^ 

exp(#) = exp(-i J B(P x )(T 1 - T 2 ) + zP x (Xi - X 2 )) (2-70) 
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= exp(im^/(Ti - T 2 ) 2 - (A^ - X 2 ) 2 ) 

2 |7\ -r 2 | , 

= exmim ) 

E(P X ) 

In the above equations, a transverse component and a longitudinal component of the 
momentum are defined as Ea. (|2-43|) . The variances in the transverse directions and 
the longitudinal direction, oj< and ctl, are given as 

GT = 1 + ( ^ )2 ( 2 - ?1 ) 

<7L= 1 + f m-2b)m^ 2 - (272) 
V 2 ff E(P x ) a ^ 

Thus the variances or and <7£ decrease with time. The spatial extension of the wave 
packet in the transverse direction is inversely proportional to err 1 / 2 and the extension 
of the wave packet in the longitudinal direction is inversely proportional to aj}l 2 . 
Hence the wave packet expands differently in two directions and the expansion in the 
longitudinal direction is determined by the absolute value of mass. For the massless 
case or very small mass case, a wave packet does not expand in the longitudinal 
direction, since the velocity is the constant in this case and ai is equal to a. 

The stationary momentum Px is proportional to the direction X\ — X 2 . After 
the wave packet expands with time, the wave becomes almost spherical wave which 
is a linear combination of many plane waves. When the measurement is made at 
the position X2, the corresponding wave has a wave vector which is proportional to 
X\ — X2 and the phase which is proportional to the proper time in the asymptotic 
region. 

In the asymptotic region, \T\ — T 2 \ — > 00, the matrix element behaves as, 
d z p{P l ,X l \^e i (p\P 2 ,X 2 ) = Nexp(i<f>), (2-73) 

* - &^W^ rin -*-§«* - A)2 + A - A,2,) ' 

(2-74) 



exp(i^) = exp(imy (Ti - T 2 ) 2 - (Xi - X 2 ) 2 ). (2-75) 

The probability to find a particle with a momentum P2 at a position X2 is given as 

P = dP 2 dX2j^\N\ 2 exp(-a{(P x - A) 2 + {Px ~ P2) 2 }). (2-76) 

For a massless case or an extremely small mass case, in large time of \T\ — T 2 \, 
the matrix element Ea.( l2-73*J) is replaced with, 



/ 



„ E(p)(T 1 -T 2 ) -> -. 

d^p{P x ,X l \p)e 1 (p\P2,X 2 ) = iVexp(#), (2-77) 

,i(T x -T 2 ) x a 5 5,2 , /5 n\2i 



N = Q U ) exp(--{(P x - P,y + (Px - P2) 2 }), 
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§3. Generalized scattering matrix 

We study many particle systems where one particle wave functions are described 
by wave packets. A field operator is expanded by a complete set of wave packets and 
coefficients become operators. 

3.1. Expansion of field 

A field operator 4>{x) is expanded by a set of momentum eigenstates as 

<j>(x,t) = J dp({p\x)-^=a(^,t) + {x\p)-jL=a\p,t)). (3-1) 

Operators a(p,t) and its conjugate a(p, ty are creation and annihilation operator of 
the momentum p. 

The operator A(Pq, Xq, Tq, t) and its conjugate are defined as linear combinations 
of the operators a(p, t) and its conjugate, 

A(P ,X ,T ,t) = J dpa(p,t){p\P ,Xo,T ) (3-2) 
A\P ,X ,T ,t) = J dp(P ,X ,T \p)a\p,t). (3-3) 

Conversely a(p, t) and its conjugate are solved from the above equations by using 
Eq. flgH ) as 

a(p,t) = J ^^A(P ,X ,T ,t)(P ,X ,T \p) (3-4) 
a\p,t) = J ^^(p\P ,X ,To)AHP ,X ,T ,t) (3-5) 

The wave packet size a in in-state is determined from the beam and the a in out-state 
is determined from the detector and they are different generally. For simplicity, we 
use the same value and omit to write in most parts of the present paper. In several 
places we write a explicitly. 

By substituting the above expansions to Eq. (j3-2jl . we have 



/dP' dX 1 -> — - -> - — 

(27r) 3 ^0' ^ol-foi Xo, To) 

(3-6) 

/rJP'd X' 
^3S(Po ! Xo,r |^,X^r^t(P ',X^^,t). 

(3-7) 

These equations show that in the space of operator A(Pq, Xq, Tq, t) the transforma- 
tion function (Pq, Xq, Tq\Pq, Xq, Tq) plays a role of Dirac's delta function. 
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By substituting the above expansion to Eq.JEEJ), the same field operator is ex- 
panded by a set of minimum wave packets of the center momentum, Pq, the center 
coordinates, Xq, and the time, To as, 

<K2, t) = J d ^*° (C(P ,X , T ;x)A(P ,X ,T , t) (3-8) 
+C(P ,X ,T ;x)*A^P ,X ,T ,t)), 

C(P ,X ,T ;x)= [ - 7 fi=(P ,X ,T \p)(p\x) (3-9) 

is obtained. It is convenient to define transformation matrices, 

C(P , X ,T ; x) = J dp^/M£j(P , X ,T \p)(p\x). (3-10) 
These matrices satisfy 

/HP rl^C 
1 (P , X , To) A, Xi, T 1 )C{P 1 ,X 1 , Ti; 3) 

(3-11) 

/dP dX 
1 1 (P , X , T \Pi, Xi, Ti)C(A, X U Ti; S) 

(3-12) 

y dxC(P ,Xo,T ;x)(C(P ,X ,T ;x)r = (Pq, X , T \Pi, Xi,Ti). 

(3-13) 

3.2. Complete set of many body states 

Many body state is constructed by the second quantized operators and the vac- 
uum. In the momentum representation the operators a(p) and a^(p) satisfy 

[a(p,t),a^,t')}6(t-t') = 5{p-$)5{t - t') (3-14) 
a(p,t)\0) = 0. (3-15) 

A complete set of many body states are constructed from 

|0), \pi), \pi,P2), \px,P2,Pz), -, \Pl,P2, -, ~Pn), (3T6) 

where 

\p 1 ) = a\p u t\0), (3-17) 

\P1,P 2 ) =a t (pi,t)a t (# 2 ,t)|0), (3-18) 

\pi,P 2 ,-,-PN) =n l at(p l ,t)\0). (3-19) 

These particle states satisfy orthonormality conditions with Dirac delta function, 

(Pi\Pi) =m~Pi), (3-20) 
(Pi)P2|Pi,P2) = &(Pl -Pi)$(P2 —PV) + permutation. (3-21) 
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Let define a projection operator / as 
I = 

|0)(0|+ J dpa\p)\0)(0\a(p) + j ^^o t (pi)at(p' 2 )|0)(0|a(pi)a(p 2 ) 

+ J -^n l dpia f (p l )\0)(Q\n l a(p l ) + • • • , (3-22) 
and multiply 7 to a state 

m = J fifing (fi)\Q). (3-23) 
Then from Eq. (|3-14|) and (|3-15j) . we have 

i\v) = J n^qrfiqi) j j^dp^mWim -pj) 
= j fmn^mo) 

= \&) (3-24) 
Hence the completeness condition, 

I = 1 (3-25) 

is satisfied. The time t is arbitrary and is omitted in the above equations. 
Similarly the operators in the mixed representation satisfy 

[A(P , Z ,T ,t),A\P , X' X, t')W ~ = <^o, X , T \P Q , X' , T )6(t - t') 

(3-26) 

A(P o ,X o ,T ,t)\0) =0. (3-27) 

The particle states defined by these operators are normalized and are not orthogonal 
even though momenta and positions are different. However the same complete set is 
constructed by the vacuum and creation operators in the mixed representation. Let 
define an projection operator I', 

/' = 

|0)(0| + J ^^aHP,x,t)\o)(o\a(P,x,t) 

+ j^n l ^^n l A\P l ,X h T l )\0)(0\n l A(P l ,X l ,T l ) + --- (3-28) 

The time t in the operators is arbitrary and is omitted in the above equations. Let 
multiply the operator I' to an state, 

\*') = j n^^^F{Qi,%Si)niA\Q h Y h Si)\Q). (3-29) 
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Then from Eq. (j3-26|) and (j3-27j) . we have 

= / n l d ^F{Q i ,f i ,s i ) J ^n m ^^A\P m ,X m ,T m )\0) 
(0\A(P m ,X m ,T m )niA\Q h Yi,Si)\0} 

= J nt^^FiQ^SMHQuYuT^O) (3-30) 

= W)- 

Hence the completeness condition, 

I' = 1 (3-31) 
is satisfied. The time t is arbitrary and is omitted in the above equations. 

3.3. Time evolution 

A unitary operator which translates a time of field operators by a finite value, 
t, is given by a Hamiltonian H as, 

U(t) = e M r. (3-32) 
In a free scalar theory the Hamiltonian is given as, 

H = J d 3 x(i(vr(x) 2 + \{V<t>{x)f + \m 2 <p(x) 2 ). (3-33) 

A commutation relation between the field operator and its conjugate, 

\<j>{x), ir(y)]5(x - y ) = i5^(x - y) (3-34) 

leads the operators a(p) and o){p) in Eq. ([3-l|) satisfy the equal time commutation 
relation Ea. (|3T4|) . The Hamiltonian is expressed as 

H = J d 3 pE(p)(a^p)a(p) + 1) (3-35) 
E(p) = \f p 1 + m 2 

and satisfy 

[H,a)(p)}=E{p)a)(p) (3-36) 

[H,a(p)] = -E(p)a^p) (3-37) 

From the above commutation relations, we have the time dependence of the creation 
and annihilation operators in the momentum space, 
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4. 4. E(ff)t , 

= U(t)a\p,0)U\t) = e^a\p,0) (3-38) 
a(jp, t) 

= U(t)a(p,0)U\t) = e ^a(p,0). (3-39) 

The creation and annihilation operators in the momentum space change only the 
c-number phase with time. The states created by the creation operators in the 
momentum space stay in the same state, 

U(t)\&) =e^ z -r J -*|Sr) (3-40) 
\9) = tJ(fr,0)<f(fr,0) a\p h 0)\0). 

The creation operators and annihilation operators in the mixed space satisfy 
a commutation relation Ea. ()3-26(l . The overlap function in the right hand side is 
computed in the next subsection. This function does not vanish even though the 
momentum or the coordinates are different. Hence the operator of one set of center 
coordinate and center momentum do not commute with the operator of a different 
coordinate and a different momentum, but a set of operators satisfy completeness 
condition from Ea. (|3-31|) . The operators in the mixed space evolve with time as, 

A^(P ,X ,T ,t) (3-41) 

= U(t)A\Po,X ,T ,o)uHt) = J d z pe^{P Q ,X Q ,T G \p)<J(p,Q) 

J d 3 pe^(P J ,To\^(^An}AHP^X'o,T ,0) 
A(P ,X ,T ,t) (3-42) 

- -> j. f o E(p)t 

U(t)A(P ,X ,T ,0)uHt) = J d 3 pe-^(p\P ,Xo,T )a(p,t) 

j d 3 pe-^(P^JiT \p^(m,X ,T )A(P^J' ,T ,0)^ 

The time dependent phase factor is not factored out and the states created by the 
operators in the mixes space, A^ (Pq, Xq) change with time, 

U(t)\V') ± e^-H-*^) (3-43) 
\V) = A\P l ,X 1 )A^P 2 ,X 2 ) At(P z ,^)|0). 

The matrix elements in the above equations are obtained in the next subsection and 
we will see that these states are approximate eigenstates. 

3.4. Generalized scattering amplitude and transition probability 

A scattering amplitude where particles in the initial state of momentum P£ are 
prepared at positions X\ and times Tf and particles in the final state of momentum 
P£t are measured at positions X^ and times are described as, 

Sout,in (3-44) 



18 



K.Ishikawa and T.Shimomura 



= (PrX^T ,; -; P L X° L T° L] out\in; P\X\T\- -; I» M X X M T* M ), 

\in- PiXlTr, -; 44f M ) (3-45) 
= A\PlX\Xi)A\Plxin) - -A\Pi,X},Ti)\0) 

\out-P^XlT 1 ---P M X M T M ) (3-46) 
= tf{PZ,Xl,TZ)A\P$,XZ,TZ) ~ -A\P?,X?,1f)\Q), 

T{ <C T t ° (3-47) 

The differential transition probability from an initial state to a final state is given 
as, 

dP = nf~i f * |^ottt,m| 2 ) (3-48) 

(Z7TJ 

and the total transition probability is obtained by integrating positions and momenta 
as, 

1 l=L dPidXi 2 

~L\ 1=1 (27r) 3 ou * ,m ' (3-4yj 
In the interaction picture, the S-matrix element is computed by 

mt^MPuXuT^Te^ J dt'^^n^ =1 A\P m ,X m ,T m )\0). (3-50) 

As will be seen in examples of the next section, the differential probability de- 
pends upon the sizes of wave packets in the initial states and final states. But total 
probabilities become universal values that are independent from the sizes of wave 
packets owing to the completeness of the many body states. 

Firstly, the total probability from one initial state to a final state of fixed number 
of particles become independent from the wave packet sizes of final states. To see 
this, let define a generalized S-matrix of wave packet size a Q of the final state and Oi 
of the initial state, 

= {P^X{T\, a - -■ PlX° L T° L a ; out\in; P\X\T\a % ; -; P^X^T M °i) ■ 

Using complexness relation of wave packets for an arbitrary wave packet size , 
the total probability from one initial state to a L-particle state of one value of a Q 
satisfy 

PM(<T ,<Ti) (3-52) 
= Jj n i=o J < ^^(in]ai\S\L,P l ,Xi,out]a )(L,P h Xiout]a \S^in-,a i ) 
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Thus the total probability is independent from the size o~ . The above total proba- 
bility also agrees with the total probability from the initial state of wave packets to 
the final states of momentum eigenstates, 



PmOo,^) (3-53) 
= Yi n l=o J j^^j^( in ^i\S\L,pi,out;)(L,pi,out]\S^\in;ai) 
= PM(momentumstate, a{). 

In the above derivations, we have used the fact that the set of N-particle states 
is complete regardless of the wave packet sizes. Thus the scattering amplitude for 
any complete set of functions gives the same total probability. Conversely we can 
compute the total probability by using the momentum eigenstates for the unobserved 
particles as far as the boundary conditions are satisfied. 

Secondly, the total probability from one state to any possible final states, 

^P(m^L) (3-54) 



£ J7 n '- 



LI <=u J (2tt) 3 
in) 

1 



U ' IXl (in\S\L, P h X h a ; out)(L, P h Xi,a ; out\S^\in) 



becomes unity. The facts that the particle states are normalized and S is unitary 
are used in the above derivation. Thus the standard probability interpretation for 
the square of the absolute value of the amplitudes is applicable with the phase space 
defined in Eq. ()3-55|) despite of the non-orthogonality of the states. 

Finally an inclusive probability where a partial set of kinematical variables are 
measured and other variables are unmeasured satisfies also a similar universal rela- 
tion as the total probability. Namely this probability depends on the size of wave 
packet of measured particles but does not depend upon the sizes of wave packets of 
unmeasured particles, if the different values of wave packet sizes are used. 

It is summarized as follows: The probability depends upon the sizes of wave 
packets of measured particles and the probability does not depend on the variables 
of unmeasured particles, such as the momenta, positions and wave packet sizes. 



§4. Few body scattering 



Few body scattering are studied as examples. We study a scattering process 
where a particle of a momentum Pi which is prepared at a space and time coordinate 
X\ , T\ and another particle of a momentum P2 which is prepared at X2 , T2 collide 
and a particle of a momentum P3 at X3, T3 and another particle of a momentum P4 
at X^T± are measured first. 
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A system with an interaction Hamiltonian, 

H int = f dxj^x) 4 (4-1) 

is studied. This interaction is of short range, and the generalized amplitude shows 
characteristic dependences on the momentum as well as coordinate. 

4.1. effective sizes of the interaction region 

Let substitute the expansion of the field Eq. (|3-8|) and the interaction Hamiltonian 
Ea. (|4-l(l into Eq. ()3-50j) . Then we have the scattering matrix in the first order of Hi n t, 

mf^MPiiXuTt) J dt '^^ln 2 m=1 A\P m ,X m ,T m )\0) 

= \J dtdmC(Pi,X u T u x,t)n m {C(P mi X m ,T m ,x,t)y. (4-2) 

We have used the two point function of the field of mixed representation and the 
field of coordinate representation, 

(0\4>(x,t)AHP,X,T)\0) 

= J ^C(Po,Xo 1 x,i)(0|i(Po,io ) t)it(P ) X,T)|0). (4-3) 

This is written further by combining Eq. (|3-8j) and Eq. (|3-ll|) as, 

(O\0(x,t)A\P,X,T)\O) = C(P,X,T,x,t). (4-4) 

In the parameter regions we are interested in this paper, this function and its partner 
is approximated well with a very good accuracy as, 

C(P, X, T, x, t) = jl (P, X, T\x, t) (4-5) 

2E(P) 



C(P,X,T,x,t) = yj2E(P){P,X,T\x,t) (4-6) 

In the following calculations we use these formula. 

The coefficients C(Pi, X[,T[, x, t) and their complex conjugate give the values of 
wave functions at (x,t). In the region where times 7] are close to the t, the product 
of the functions is a Gaussian function around the peak in the variables x and t and 
is expressed as 



(N 3 ) L nexp(- — (x-X J - Vj(t - Tj)) 2 ) x 



n t c(Pi,x h T h £, t)n m (c(p m , x m , r m , x, t)f (4-7) 
1 

exp (-iE{Pi)(t - Ti) + iP l ■ (x-Xi) + iE{P m )(t - T m ) - iP m ■ (x - X m )) 

ATexp(- - f ) 2 - J-(t - t ) 2 y*. 

10~s Ao~t 
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The peak position xq and to are determined as, 



xo = a s y^—xi(t), 

^ CFl 

B 



(4-8) 



(4-9) 




x l (t)=X l +v l (t-T l ) 




(4-11) 



and the variances as and ot are determined as 

a 



(4-12) 



-l 



(4-13) 



l i 



The region around the peak within the spatial distance \/a s and the time distance 
\fa T gives a dominant contribution to the integral. Both distances are proportional 
to the a. So if the as is small, the ax is also small. The normalization N and the 
phase 4> are complicated functions of the energies, momenta, spatial positions, and 
temporal positions. Explicit formulas are given in Appendix. 

4.1.1. Complete measurements 

So far, all particles have the same wave packet size and all particles are measured. 
When this is not hold and different particles have different wave packet sizes and 
some particles are not measured, the effective sizes as and ax become different from 
the above values. We study the behavior of these variances in general cases where 
each wave packet has its own size and all momenta and positions are measured here. 
Let specify the wave packet size of the 1-th particle as a; , and its velocity as vi , then 
the variances as and ax are given by 



The as is determined mainly by the small 07 of the measured particles but the ax 
is determined by the a% and v% of the measured particles. The small a\ does not 
contribute if the corresponding v% vanishes. So the time variance ax could become 
large even though the space variance as is small. The large a t is compatible with 
the small as in this case. 

The vi depends on the momentum Pi. So cjt is not a constant generally but 
varies in the kinematical region of the final state. The exception is the case when 07 
is infinity, i.e., plane wave. 




(4-14) 



(4-15) 
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4.1.2. Partial measurements 

When some portion of particles are measured and others are unmeasured, the 
probability depends on the wave packet sizes of the measured particles. When one 
particle, I = 1, is measured and other particles are unmeasured,it depends on the 
wave packet ai. Hence the effective sizes of the vetcex area for computing the total 
probability are obtained by letting 07 = 00 for I 7^ 1 

a s = a 1 (4-16) 

a T = (—v 2 - —vfr 1 = 00 (4-17) 

The spatial size a$ is determined from the o\ of the observed particle but the tem- 
poral size ot diverges. 

Next, the probability when two particles are measured and other particles are 
unmeasured depends on the wave packet o\ and 02 and velocities v\ and V2 of the 
measured particles. The effective sizes are obtained by letting 07 = 00 for Z 7^ 1, 2 



a s 



ax 



E-M =^r (418) 



a\ + a 2 
1=2 



yZ-vf-asiyZ-vlf) = 7 ^^ (4-19) 



\l=l 

aT diverges if v\ is equal to V2 ■ 
4.2. Short distance scattering 

When two particles are in the initial states and two particles are in the final 
states and the distances |JQ — XA and |Tj — Tj\ are small, the formula Ea. (|2-28|) 
for the small time differences is used. We expect that the amplitude in this region 
shows features of the translational motion of wave packets and other features of the 
generalized scattering amplitude. For simplicity, we present the results when all the 
particles have the same wave packet size. The general case is given in the Appendix. 

The transition amplitude in the lowest order of Hj n t is given by 

A J dtd^ l C(P l ,X l ,T l ,x,t)n m (C(x,t,P m ,X m ,T m ,x,t))* (4-20) 

= XNln^EiPt))- 1 ' 2 J dtdxe* V Y,{-isi{{t ~ T{)E{Pi) - P t • (x - X)} 

ex P -X-m- Ti)} 2 }, 

= Nexp(i<j) - R), 

N = ( - ^ _ 2 ) 3 / 2 77;(2g(P,)r 1/2 , (4-21) 



sgniiTiEiPi) - Pi ■ X { ) + 4>, (4-22) 



Generalized S-matrix in Mixed Representation 



23 



R = |(E + ^ - ^ - + 

where sgrtj is a signature 

s^nj = +l(in — state), —l(out — state) (4-24) 

and and R are small quantities and are given in the appendix for most general 
case. In the above equations, the dominant part in the phase, <f>, is the standard 
phase of the plane wave of the momentum p*o and energy Eq. Prom the first factor 
and third factor of normalization R, the momentum conservation is approximately 
satisfied with the variance, (f ) -1 / 2 and the energy conservation is also approximately 
satisfied with the variance — ?£- ^ y > )~ x l 2 . The second factor of the R shows that 

the particle trajectories coincide and coordinates X{ — v(Ti are the same within the 
distance (8a) 1 / 2 . So particles follow classical trajectories. 

From the amplitudes, we define new transition probabilities that depend upon 
the coordinates in addition to the momenta and argue on an asymptotic condition 
of the standard scattering amplitude which depends upon the momenta. 

The transition probability is a square of the absolute value of the amplitude and 
is expressed as 

P(P 3 ,X 3 ,T 3 ;P A ,X 4 ,T A ) (4-25) 
1 1 



(JV) 2 exp(-^p. )2 



2!(2tt) 3V ' ^ 4 

This has a peak at the positions where the conditions 

Xi - - Xj + VjTj = (4-26) 

are satisfied. Thus the peak is along a line 

Xi = v^ + C, (4-27) 

where C is a constant vector. These positions depend upon the times and the time 
Tj are arbitrary in the present formalism, hence it is possible to choose the times in 
such manner that these positions are inside of detectors if the detector is located in 
the direction of the momentum. Then total probability integrated on this direction 
is measured. To see this probability, let us decompose the position vector of the 
i-th particle into the longitudinal component and the transverse components with 
respects to the velocity v 1 , 

X 1 = v l S i + 4lV, (4-28) 
v 4 ■ r? T = 0, (4-29) 
n % T -n! T = bij. (4-30) 
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The volume element is written as 

dX j =dS s dX i T \v i \. (4-31) 
Using these variables, the Gaussian factor of the differential probability is written 

as 

-tfVP -X j +v j T j f (4-32) 

u 

= ^(^(S* " T ' 1 ) + ™tXt ~ &(S j ~ T j ) - 4^t) 2 
u 

= J^S* + rf T Xf - W - 4^t) 2 > 
u 

S i = S i - T\ (4-33) 

Thus the longitudinal variable is combined with time T % and it is possible to 
replace the longitudinal coordinate with the time variable in Eg. (J4- 34(1 , 

The transition probability is given by an integration of a differential probability 
over the momenta and coordinates as, 

P = J n™=I^^P(P 3 ,X 3 ,T 3] P 4 ,X 4 ,n). (4-34) 

In the ordinary detectors neither the the precise value of the time T l nor the lon- 
gitudinal coordinates X 1 l are measurable but the total probability P(P 3 ,Xj,, ; P 4 , Xj) 
integrated on these variables is measured. 

To obtain this probabilitythe variable Si is integrated. Then the probability 
P{P 3 , Xf,, T 3 ; P 4 , X|, T 4 ) is found and is written as 

P(P 3 ,X^,T 3 ;P 4 ,X^,T 4 ) (4-35) 

= J 7i™=3 4 dxr^^p(p 3 ,x 3 ,r 3 ;j p 4 ,x 4 ,r 4 ) 

= P(P 3 ,X 3 ,;P 4 ,X 4 ). 

The time dependence disappears in P(P 3 , X|,, T 3 ; P 4 , Xj, T 4 ). 

Next, we make a connection of the present result with a standard scattering ma- 
trix where an asymptotic condition is satisfied and only the momenta are observed. 
In the ordinary scattering processes the initial time Tj is — oo and the final time is 
+oo. The distance \X,i — Xj\ in the initial state is proportional to \vi — Vj\Ti and the 
distance |Xj — Xj\ in the final state is proportional to \vi — Vj\Ti. They become large 
except \vi — Vj\ = 0. When we define this case from a limit \vi — Vj\ — * where a 
large T limit is taken first, a distance between two wave packets becomes infinity and 
the wave packets at T — > ±oo do not overlap each others. The theory thus defined 
satisfies asymptotic condition. 
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By integrating the coordinates, we have the momentum dependent differential 
probability, 

P(P 3 ,;P4) = Nexp(-jC£P t f - t a (^i?^) 2 ) (4-36) 

4 ^ i z2i,j( v i- v j) ^ 

N = l(4a7r) 3 2- 3/2 exp(- — (Xi - wiTi - X 2 + 77 2 T 2 ) 2 ). (4-37) 

Z 0(7 

In Ea. (l4-36l) . the normalization factor N is a constant which does not depend on 
the final state and the momentum dependent probability is almost the same as the 
probability of the plane waves. By integrating momenta of the final states,we have 
the total probability. 

If the initial state and the final state have different values of the g, we use the 
formula given in the appendix. Let o~ be the size for all final particles and the Gi be 
the size for all initial particles, the probability is given as, 

P(P 3 ,;P4) = iVexp(-a s (^. J P i ) 2 -a t (^^^) 2 ) (4-38) 
N = ^(4vr) 3 (^) 3 / 2 exp(--^-(X 1 - - X 2 + v 2 T 2 ) 2 ), (4-39) 

°s = \(- + -)-\ (4-40) 

a t = (ufi-%-\ (4-41) 

G 4 G. q 



4.3. Long distance scatterings: first order 

When one of the times, T\, is in a position far away from other times T}(Z ^ 1), 
and the classical trajectories meet at around a time near Ti(l ^ 1) and coordinate 
X[(l ^ 1), the dominant contribution in the integration comes from the region near 
Xi(l 7^ 1) and the position Xi — v{Fi. One of the time difference t — T\ becomes 
large and asymptotic expansion for the corresponding C(Pi, X%, Ti\t, x) is used. The 
transition matrix element becomes, then, 

mf^AiPuXuTt) J dt'^^nl =1 AHP m ,X m ,T m )\0) 

= A ( =, Jl =, f (^JV^e^^-^^ 

2E(P 1 )2E(P 2 )2E(P 3 )2E(P i ) J y 

jy 3 * e -^^-^2-v2(t-T 2 ))HiE(P 2 )(t~T2)-iP2(x^X2) ( 4 . 42 ) 
x n j=3A N 3 e- h (*~*J (*- T ^) 2 -iE(P 3 )(t-T 3 )+Pj (x-Xj) 

where N asym and the stationary momentum P-g i _ s are 

Nasym = M^-^ (^-^(^ (4-43) 
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Xl ~* "'((Ti -tf- {X x -x) 2 ) 1 / 2 

Substituting these expressions, we have 

= A( — ^— -5 L^TT^ ) 1/2 (l^3| 2 ) 2 (^— ) 1/2 ( ' 



2E(P 1 )2E(P 2 )2E(P 3 )2E(P 4 )' " '■'^ 2 hk + l J ^2^r +1 ' 

v 2 (t-Ti) 2V ; 2 m 4 (i-Ti) 2V ; ; 

tfe* exp (-^L(x - x' ) 2 - ±-(t - t' ) 2 ), (4-44) 
where the center position xq is given by 

x Q = X l + (t-T 1 )-^- (4-45) 

and the normalization , the phase, the variances, and center positions are defined by 
diagonalizing the products of wave packets, 

N e i exp (- — (£- 4) 2 - —(t - g 2 ) = N3 * e -^-^-Mt-T2)r+iE(^)(t-T 2 )-i^-X 2 ) 

2a s 2a t 

x/T^^e-^^-^-^^-^W'^^^^^^-^J. (4-46) 
If the variances <r s and at are small values, the amplitude is further written as, 

w 1 ) l / 2 (\N \ 2 ) 2 ( 1 ) 1/2 ( 1 ) 

K 2E{P l )2E{P 2 )2E{Pz)2E{P i ) > U 3|jl2 m + i j ^ + 

iVe^(2o- s 7r)l(2at)i (4-47) 

This expression of the amplitude shows that the wave packet expands and has the 
phase factor which is proportional to the square root of the proper time. 
When we integrate on the variables x, t first in Ea. (|4-42|) . we have, 

A J dtdxn l C(P h X h T l \x,t)n m C(P m ,X m ,T m \x,t)* (4-48) 

= (|iv 3 | 2 ) 2 &/ 2 (— ^--^) 1/2 ^ 
3 < [vy > — < v > z 

I E ( 449 ) 
l E ^-) 2 ' ( 4 - 5 °) 



6 ,=1,3 



< -i; 2 > 



3 ^1,3 
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where 5Ti, SXi, and 8P\ are of the order 0(1). N is a normalization factor which 
depends upon kinematical variables. The integration on the variable p is carried 
with a use of the stationary phase approximation as in the previous cases. 
Finally we have the amplitude, 



a J dtdxniC(Pi,Xi,T l \x,t)n m c(P m ,x m ,T m \x,ty 



Thus the amplitude depends on the large variables X\ and T\ in a simple form. 
The normalization factor is inversely proportional to T\ and the phase factor is 



proportional to the mass and the proper time my c 2 (Ti) 2 — [X\) 2 . 

4.4. Long distance scattering :second order 

Next we study the few body scattering amplitude in the second order of interac- 
tion where there is one propagator D(t,x). The propagator connects two interaction 
points, (xi, ti) and (x2,t2), 

mf^PuXuTt) J dt 1 dfa r( ^ ( ^ g ^ (fa)) /z^ 1 At(P ro ,x ro ,r ro )|o) 

= A 2 y dtidxi J dt 2 dx 2 V g (t 1 ,x 1 ,Pi,- ■ ■)D(ti-t 2 ,x 1 - x 2 )Vg(t2,x 2 ,P2,- ■ ■)*, 

V g {t 1 ,x 1 ,P h ---) = n l C{P h X u T l \x 1 M), (4-52) 
V g (t 2 , x 2 ,P 2 , ■■■)* = n m C(P m , X m , T m \x2,t 2 )* . 

The propagator D(t± — t2, x\ — X2) is given by, 

D( tl - t 2A - ft) = i I J ^^M'>-*X m =^- (4'53) 

We study the configuration when times = 1,3) are close each others and times 
T m (m = 4,6) are close each others but the first group of times is separated from 
the second group of times with a large distance. Regions when the time variable 
t\ is near Ti(l = 1,3) and the other time variable t 2 is near T m (m = 4,6) or the 
opposite give the dominant contribution to the amplitude in the time integration. 
Since the distance \ti — £2) is large, on mass shell kinematical region where p 2 = m 2 
is satisfied is dominant in the momentum integration. Using the stationary phase 
approximation in the momentum integration, we replace the propagator with the 
asymptotic form that is obtained at the stationary momentum p x 

Dasymih - t 2 ,X 1 - X 2 ) = iN x -^±^ j^-**\ E{j(x) (4-54) 

N x = (— ) 1/2 (— ), (4-55) 
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1 (h ~ h) m 2 

^ = 2-m^ nL=7T m^> (4 ' 56) 



where the momentum p x is given as, 

m 

p s ={x 1 -x 2 ) — — — 4-57 

((*l - h) 1 ~ (xi - x 2 ) 2 ) 1 / 2 

We substitute these expressions into the amplitude and we have, 

X 2 j dtidxidt 2 dx 2 V g (ti, xi, Pi,)D(ti -t 2 ,xi - x 2 )V g (t 2 , x 2 , P m , X m , T m )* 

= \J dhdxxVgit^xx,^ X u T iY PxXl \E{p*)iN x { \ E( ^ ) 1 ' 2 (4-58) 

A j dt 2 dx 2 V g (t 2 , S 2 ,P m ,X m , r m )*e~^ | E(fe) ( ,^1^ f' 2 - 

In integrating (ti,Xi),i = 1,2 the momentum p x of Ea. (|4-57[l becomes a constant 
vector, if these variables are in the narrow regions. Then the total amplitude is 
proportional to the product of the two amplitudes and is inversely proportional 
to the spreading of the wave packet jt ■ The total probability is proportional to 
the product of probabilities of two processes and is inversely proportional to 
Interference effect of propagating wave is negligible in this case. Thus in the present 
regime the intermediate state is treated as an observed particle. Hence the single 
particle treatment of the intermediate state is applicable. Using the argument of 
Ref. 7 ) and, 6 ) the effect of the wave packets are described by the ensemble of the 
energy eigenstates in this regime. 

In an opposite situation where these integration variables cover wide regions 
and p x is not a constant vector but varies with these variables (ti,x\) or (t 2 ,x 2 ), 
the naive single particle treatment is not justified. The total amplitude becomes 
linear combinations of the amplitudes of the various values of the momentum, p x , 
and become different from the product of two amplitudes. We split the integration 
regions into small regions V\ and obtain the momentum defined from a pair of 
these regions. Using them we have the amplitude, 

X 2 J dtidxidt 2 dx 2 V g (ti,xi,P h )D(ti - t 2 ,xi - x 2 )V g (t 2 , x 2 , P m , X m , T m )* 

= E A j Vi dhdx x V g {t u x u P h Xl, T|)e^ ( 27r )32^/) ^ (4 ' 59) 

A J dt 2 dx 2 v g {t 2 , x 2 ,P m , x m , r m )* e-*^ \ m , ) ( {2lTf l m) ) 1/2 - 

The total amplitude is a linear combination of the amplitudes of different momentum 
$ x . Probability may show the interference of the different intermediate momentum. 

In this situation we are able to write the amplitude in a different manner. Let 
write the propagator as, 

D(ti -t 2 ,xi - x 2 ) = -2i J d 3 xD(t 1 -t,x 1 - x)D(t -t 2 ,x - x 2 ). (4-60) 
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We substitute this expression and we have the amplitude, 

X 2 J dhdx^dxiVgih,^, Pi,)D(h - t 2 ,xi - x 2 )V g (t 2 , X 2 , Pm i X rn , T m ) 

= -2i\ 2 j dx j dt 1 dx 1 V g (t 1 ,x 1 ,Pi,X l ,T l )D(t 1 -t,x 1 - x) (4-61) 

x j dt 2 dx 2 D(t - t 2 ,x- x 2 )V g (t 2 ,x 2 ,P m ,X m ,T m )*. 

Using the stationary phase approximation in the momentum integration, we 
have the propagators, 

Dasymih -t,X 1 -x)= iNW 1 ^\^)\ (4 . 62) 

(2ir) 3 2E(p> x >) {Px ' 

= (4y ) i/2 (4y )» ^ 



and 



(1) l (ti-t) (!) (i) m 2 

= (xi - x)— — — ^ — — - (4-65) 

((ti -t) 2 - (xi -x) 2 ) 1 ' 2 



D asym {t - t 2 ,x-x 2 ) = z 2 AT( 2 )^_L_ e ^ 2) (— 2)| Btfg0) (4-66) 

^ = (4) ) 1/2 (-^) )• ( 4 ' 67 ) 
«7i !7t 



(2) _ 1 (t - tg) (2) _ (2) rn 2 



3° = <'-*> ««-w» -?*-*»■/» • < M9 > 

The amplitude becomes, 

-2iA 2 y dxT^XuPu ■ ■ ■ , ; x, t)N^N^T 2 (X m , P m , *)> 



(4-70) 

T^XuPu--- ,;f,t) = ^-3 (4-71) 



1 

2(2^) ; 



x f dhd^VgituX^Vji^—tj—e^^-^l 



1 

2(2^) ; 



T 2 {X m ,P m ,--- ,;x,t)* = 3 (4-72) 
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This amplitude agrees with the previous form if the momentum j4 and pi are 

regarded as constant vectors. If these momenta are not constant vectors, the x 

— -ii) 
dependence and Xj,j = 1,2 dependence of p£ ,j = 1,2 are taken explicitly in the 

time and coordinate integrations. 
4.5. Factorization 

Factorization is a general feature of the amplitudes obtained in the previous 
sections. Namely the amplitudes are factorized into amplitudes of sub-processes 
that depend on close space and time coordinates X^jT^ where — Tj x \ « 
and Xi 2 ,Ti 2 where |Tj 2 — Tj 2 \ and \Ti x — Tj 2 | oo. In these short distance 
amplitudes, Gaussian momentum integration around minimal are applied and am- 
plitudes become almost equivalent to ordinary scattering amplitudes. On the other 
hand, the long distance parts have particular forms that are proportional to the 
inverses of time difference \T^ — Tj 2 | in the small mass case and to the phases 

exp(im^ (T{ 1 — Ti 2 ) 2 — (X^ — Xi 2 ) 2 ) , where and X^ are average values of times 

and positions in a group 1 and Ti 2 and Xi 2 are average values of times and positions 
in a group 2. The former behavior is due to the expansion of the wave packets. It is 
possible to decompose amplitudes in general many body amplitudes in which space 
time coordinates are separated into many groups. Each amplitude for the process 
of close coordinates is almost equivalent to ordinary scattering amplitude and the 
amplitudes for the long distance parts have the particular normalization and the 
phase factor of the above forms. 

§5. Summary 

We have defined the generalized scattering amplitudes which have dependence 
upon particle's positions in addition to the particle's momenta. Idealistic cases where 
the positions and the momenta satisfy minimum uncertainty relations are studied 
by the use of minimum wave packets, coherent states. 

Since wave packets are linear combinations of eigenfunctions of free Hamiltonian, 
wave packets change with time. Wave packets move with a constant group velocity 
and expand. These behaviors occur since each wave of definite momentum has 
a different velocity. They reveal a particle's nature and a wave nature of wave 
packets. Expansion is slow and has been irrelevant to any observations in high energy 
experiments till recently. They are relevant in some long distance experiments and 
its effects are analyzed in the present work. We found also that the expansion speeds 
satisfy new uncertainty relations expressed in Ea. (|2-55|) and Ea. H2-56|) . 

Several relations which must be satisfied for the transition amplitudes and prob- 
abilities are proven. Completeness of the mixed representation is proven and is used 
for defining the weight of phase space integral for both variables of momenta and 
coordinates. The particle states which are specified by momenta and positions are 
normalized to unity and Dirac delta function is unnecessary for the normalization of 
states in mixed representation since the states are normalized but are not extended 
in space. The whole transition probability from one state to states of a fixed particle 
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number becomes an intrinsic value which is independent from the wave packet sizes 
and the total transition probability from one state to all possible states becomes 
unity under the use of the present measure of phase space even though states of 
different momenta and positions are nonorthogonal. So probability interpretation 
holds. 

Several examples in few body scatterings are analyzed and amplitudes are ex- 
plicitly computed in the lowest order and the second order of the interaction Hamil- 
tonian. Translational motions and expansions of wave packets are taken into account 
explicitly and their effects are seen in the manifest manner. It is shown that scatter- 
ing amplitudes which have long distance part in addition to short distance part are 
factorized. The asymptotic condition for the ordinary scattering, which is satisfied by 
the addition of ie in propagators in the standard s-matrix, is realized automatically 
in the present formalism by taking a suitable limit of the present amplitude. 

Applications to neutrino long distance experiments, 10 ) 13 ) and others will be given 
in separate works. 
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Appendix A 

Generalized vertices of arbitrary wave packets 

The product of the wave functions at (t, x) are the Gaussian function of the 
space time coordinates (t,x), 

n-N*e~^ ~ Vj { - t ~ T ^ )2+iE{ -^ t ~ T ^~ i ^ ^-x? } 

3 3 

x ^^ e -^(5-^-vi(t-r0) 2 -»^(p0(*-^)+»pi(5-^i) ( A .-g 

= n J N*n l N l e-^-* o{t))2 ~^ {t - to) \ R+i t 

Wave packet parameters in the spatial directions and the temporal direction are 
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^ = ^'^-^0 (A-3) 
a t (Jj J a s 

and the central values of the space-time coordinates are 

x (t) =v t + x (0), (A-4) 

v = a s £j- Vj, (A-5) 

x (0) = a s (Uj—Xj - i(Sj(±)pj)) (A-6) 
*o = 0t( — #o • £o - £j — vj ■ Xj + iSj{±)E(pj)) (A-7) 

X^Xj-vjTj. (A-8) 

The real part determines the magnitude of the amplitude and is composed of the 
trajectory terms and the energy- momentum terms. The former give constraints on 
the particle trajectories and the latter give constraints the total energy and total 
momentum and are determined as, 

R = Rtrajectory ~\~ Rmomentumi (A'9) 

1 =. 2 1 =* 2 % 2 

Rtrajectory = —Ej- Xj +2<J s (Ej- Xj) + 2a t (£j {vo — Vj ) Xj ) , 

laj zaj 

(A-10) 

Rmomentum = ~{Zj(±)(E(pj) - V Pj)f - y (27j(±)pj) 2 . (A-ll) 

The phase factor is composed of the primary term which expresses the energy mo- 
mentum dependent phase and the secondary terms which are due to finite sizes of 
wave packets and are determined as, 

= 00 + ^1, (A- 12) 

0o = Zj(±MXj ~ E(Pj)Tj), (A-13) 



1 = -2a t (^^r(^-r / )X y )fr(±)o,(P i - Edl,))) 



-2a s {Zj{±)pj)(Zj^-Xj). (A-14) 
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